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Environmental noise and disorder play a critical role in quantum particle and wave transport in complex media, 
including solid-state and biological systems. Recent work has predicted that coupling between noisy environments 
and disordered systems, in which coherent transport has been arrested due to localization effects, could actually 
enhance transport. Photonic integrated circuits are promising platforms for studying such effects, with a central 
goal being the development of large systems providing low-loss, high-fidelity control over all parameters of the 
transport problem. Here, we fully map out the role of static and dynamic disorder in quantum transport using a 
low-loss, phase-stable, nanophotonic processor consisting of a mesh of 56 generalized beamsplitters programmable 
on microsecond timescales. Over 85,600 transport experiments, we observe several distinct transport regimes, 
including environment-enhanced transport in strong, statically disordered systems. Low loss and programmability 
make this nanophotonic processor a promising platform for many-boson quantum simulation experiments. 


Quantum walks (QWs), the coherent analogy to clas¬ 
sical random walks, have emerged as a useful model for 
experimental simulations of quantum transport (QT) phe¬ 
nomena in physical systems. QWs have been implemented 
in platforms including trapped ions^’^, ultra-cold atoms^, 
bulk optics^“^ and integrated photonics^’^"^^. Integrated 
photonic implementations are particularly attractive for rel¬ 
atively large coherence lengths, high interferometric vis¬ 
ibilities, integration with single-photon sources^^’^^ and 
detectors^^, and the promise of scaling to many active and 
reconfigurable components. The role of static and dynamic 
disorder in the transport of quantum walkers has been of 
particular interest in the field of quantum simulation^^’^^. 

Control over static (time-invariant) and dynamic (time- 
varying) disorder enables studies of fundamentally interest¬ 
ing and potentially useful QT phenomena in discrete-time 
(DT) QWs. In systems with strong dynamic disorder, illus¬ 
trated in Fig. 1(a), a quantum walker evolving over T time 
steps travels a distance proportional to a/T; the coherent 
nature of the quantum walker is effectively erased, result¬ 
ing in classical, diffusive transport characteristics^^’^^. In 
contrast, a quantum walker (or coherent wave) traversing 
an ordered system travels a distance proportional to T as 
a result of coherent interference between superposition am¬ 
plitudes — a regime known as ballistic transport (see Fig. 
1(b)). Perhaps most notably, a quantum walker propagating 
through a system with strong, static disorder becomes ex¬ 
ponentially localized in space and time, inhibiting transport, 
as illustrated in Fig. 1(c). This QT phenomena is known as 
Anderson localization^"^ and it has been observed in several 
systems, including optical media^“^^’^^’^^. For systems in 
which transport has been arrested due to Anderson localiza¬ 
tion, it has recently been predicted that adding environmen¬ 
tal noise (dynamic disorder), over a finite range of strengths. 


could result in enhanced transport. This effect, known as 
environment-assisted quantum transport (ENAQT), is be¬ 
lieved to play a key role in the high transport efficiencies in 
continuous-time photosynthetic complexes^^’^^. With such 
a variety of QT phenomena accessible, it is of fundamental 
interest to develop systems capable of fully probing the dis¬ 
order space — potentially revealing new QT phenomena in 
the single- and multi-particle regimes. 

In integrated photonic systems, static and dynamic dis¬ 
order have been introduced by fabricating circuits with ran¬ 
dom parameter variations or post-processing^’^^ — render¬ 
ing explorations of large parameter spaces in QT simula¬ 
tions possible, but cumbersome. Further, single realizations 
of disorder are not suitable to characterize transport as they 
can produce a wide range of output distributions^^; ensem¬ 
ble averages over many realizations are necessary to accu¬ 
rately reproduce output statistics. While static and dynamic 
disorder have been studied separately"^’^"^^’^^, transport in a 
system capable of implementing both simultaneously — a 
requirement for the observation of ENAQT — in all com¬ 
binations and strengths, with low loss, and over a large 
number of realizations has not been demonstrated^^. De¬ 
veloping systems capable of meeting these requirements, 
and the more stringent requirement of generating arbitrary 
single-particle unitary operations'^, has been the topic of 
several recent theoretical investigations^^’^^. 

Here, we investigate the effect of all combinations of 
static and dynamic disorder strengths on the transport effi¬ 
ciency of a single-particle undergoing a DT, symmetric QW. 
Over 5,600 experiments, we probe the spatiotemporal evolu¬ 
tion of a quantum walker traversing the system under several 
distinct transport regimes. To determine whether ENAQT 
can occur in DT systems maintaining a strong static dis¬ 
order component, we perform 80,000 experiments spanning 
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FIG. 1. Conceptual representation of quantum transport experiments. Schematic of programmable nanophotonic processor. 
Horizontal lines represent nanowire waveguides; boxes represent 4-port reconfigurable beamsplitter ((see lower inset of (a); thermal 
phase modulators control the splitting ratio via ^ 1,2 and the phase shift between steps via 0i,2))- Time is defined from left to right; 
space is defined from top to bottom. Colored regions represent output differential phase setting (a) Output phases are randomized, 
resulting in incoherent, diffusive transport shown conceptually with a Galton’s board representation at the last time step, (b) Ballistic 
propagation configuration. The quantum walker is most likely found along the outer modes of the PNP in the absence of phase 
disorder, (c) Anderson localization configuration. The quantum walker is mainly confined to the central modes in the presence of 
strong static phase disorder (the phase setting on a given waveguide mode is constant in time and denoted by (j)i). 


400 unique combinations of static and dynamic disorder. 
Our data show an enhancement in the transport efficiency 
of a quantum walker in systems with strong static disorder 
as dynamic disorder is increased. 

The 85,600 experimental realizations of QWs are en¬ 
abled by a large-scale, low-loss, phase-stable programmable 
nanophotonic processor (PNP) that provides complete con¬ 
trol over static and dynamic disorder as well as coupling 
rates in QWs. Previously, photonic integrated circuits for 
quantum information processing have been limited to 30 in¬ 
dividually tunable elements and 15 interferometers^^. The 
PNP is composed of 213 individually tunable phase modula¬ 
tors and 56 interferometers spanning a chip area of 4.9 mm 
by 1.7 mm. This high component density is enabled by high 
silicon-to-silica index contrast (enabling the fabrication of 
waveguides with bend radii down to 2 jiim^^) and the large 
thermo-optic coefficient of silicon (enabling compact phase 
modulators^^). 

We consider random walks of a single photon (boson), 
which are characterized by output intensities identical to 
those obtained with coherent light. While single-particle 
statistics are easily calculated on classical computers, the 


Hilbert space occupied by a state of n identical photons 
scales exponentially with n — recent work suggests that 
even sampling the output distributions from linear unitary 
processes, including QWs, with a relatively small numbers 
of photons {n > 30 for a circuit with approximately 1000 
modes) becomes intractable for classical computers^^"^^. 
Further, QWs augmented with feed-forward controP^“^^ 
could be used in universal quantum computation schemes^^. 
These single-particle QT simulations performed using our 
large-scale, low-loss, phase-stable, and high fidelity PNP 
mark an important step towards post-classical computing 
and quantum simulation applications in the multi-particle 
regime. 


PROGRAMMABLE NANOPHOTONIC PROCESSOR 

The PNP, shown in Fig. 2(a), consists of a mesh of re- 
configurable beamsplitters (RBS). Each RBS is composed 
of two 50% directional couplers separated by two internal 
{& = 0i — O 2 ) and two external (^ = — ^ 2 ) thermo¬ 

optic phase shifters^^, as shown in Fig. 2(b, c). Redun- 
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FIG. 2 . Optical micrographs of programmable nanophotonic processor, (a) The system is 4.9 mm by 1.7 mm and is composed 
of 56 Mach-Zehnder interferometers, 213 phase shifters, and 112 directional couplers. Waveguides are outlined in black, (b) RBS 
unit cell with two directional couplers and two internal (^ 1 , 2 ) and two external ( 0 i, 2 ) low-loss, thermo-optic phase shifters, (c) 
Thermo-optic phase shifter. Heat is localized to the waveguide using narrow thermal bridges, (d) Phase versus squared bias and 
(e) transmission versus phase for all 56 internal phase shifters in the light cone. All elements exceed T.'k radians phase shift with an 
average interferommetric visibility of 0.9993 zb 0.0002. (f) Histogram of transmission minima; the majority of interferometers have 
minima at the noise floor of our experimental setup. 


dancy in internal and external phase shifters balances loss, 
increasing interferometric visibility for the single- and meta¬ 
interferometers that compose the PNP. The RBS applies 
the SU(2) rotation. 


Urbs 


sin(0) cos(0) 
cos(0) — sin(0) 


in the spatial-mode basis. Each RBS can implement any 
rotation in SU(2) by using its internal and external phase 
shifters and the output phase shifter of the preceding RBS; 
for the input RBSs in the PNP, all of SU(2) can be accessed 
by choosing an input phase. 

As shown in Fig. 2(a), the RBS lattice of the PNP is 
staggered such that a particle in the top port of an RBS 
will enter the bottom port of an RBS in the next layer. 
This nearest-neighbor coupling is described by 


.d 


. cy ^ 


where x is the spatial coordinate shown in Fig. 1(a), Cx,x+i 
iCx,x-i) is the amplitude associated with transitioning to 
a higher (lower) mode at each layer, respectively, and a is 
the mode annihilation operator. 


The thermo-optic response of all 56 internal RBS phase 
shifters within the single-input light cone used in this ex¬ 
periment is shown in Fig. 2(d): all phase shifters reach 
beyond Ztt radians of imparted phase. Transmission curves 
and histogram for the RBS unit cells are shown in Figs. 
2(e) and (f), with a mean visibility V of 0.9993 ± 0.0002 
~ i^max Imin) / {^max “b Imin)'i where I max’) I min are 

measured optical intensities in Watts). The sub-linearity in 
Fig. 2(d) is attributed to electron velocity saturation; vari¬ 
ation in phase results from contact resistance differences 
across the 240 electronic probe pads on the PNP chip. 

For bias levels corresponding to O^cj) G [0,27r], ther¬ 
mal cross-talk between phase shifters results in a transmis¬ 
sion change in nearby RBSs of less than -40 dB — the 
dynamic range of our experimental setup. The thermo¬ 
optic modulators have a 3 dB bandwidth of 130.0 ± 5.59 
kHz^^, permitting up to 10^ full PNP reconfigurations per 
second. The on-chip transmission loss through the PNP 
was measured to be less than 2.0 dB and is attributed to 
waveguide propagation loss^^. The PNP was fabricated 
in a complimentary metal-oxide semiconductor (CMOS) - 
compatible, silicon photonics process (see Methods for more 
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detail). Waveguide edge roughness results in minor differ¬ 
ences in path lengths between the arms of individual RBS 
unit cells^^. We eliminated this problem by calibrating the 
PNP and programming it to implement a known unitary 
operation (see Supplemental Information for details on cal¬ 
ibration). 

To program all 213 phase shifters, we developed a 240- 
channel electronic biasing printed circuit board (PCB) ca¬ 
pable of driving all channels from 0 to 12 volts with 183 |llV 
discretization. To measure the PNP output mode intensi¬ 
ties, we developed a 12-channel photodiode array PCB with 
approximately 40 dB of dynamic range per channel. See the 
Supplemental Information for experimental setup details. 


EXPERIMENT 


As illustrated in Fig. 1(a), we use a spatial encoding 
for the position {x) of the quantum walker with the axial 
coordinate {t) defined at each layer of RBSs. The settings 
for the QW program are shown in Fig. 3(a). In the first time 
step, we excite the bottom waveguide mode {x = 8) and 
program the first RBS as (0,^) = 1‘i) giving the 

output state (up to a global phase) |^) ^ ^(1^) + ^ l^))- 
where t (6) is the top (bottom) output mode of an RBS. We 
program all subsequent RBSs as Fladamard gates (7r/2,0), 


^ Q(^/2>7r/2) 

□ (7r/2,0),ff 


O Modulated 
0 02 Modulated 



1 1 
1 -ij ’ 

ensuring that the upward traveling terms in the superpo¬ 
sition are 7r/2 out of phase with the downward traveling 
terms and yielding a symmetric QW. 

Static and dynamic disorder are applied by randomly vary¬ 
ing ^ while holding 0 constant. Output phases are labeled 
as ^x,t, where ‘(x,t)’ is a coordinate in the space-time of 
the QW. Static disorder is controlled by selecting phases 
from the uniform distribution {U[—7r,7r]) over the position 
space of the QW ^x,t = ‘^.Tvcud * for each x, 

where cud ^ [0,1] is a coefficient weighting static disor¬ 
der strength. We consider space-time uncorrelated dynamic 
disorder which is applied as ^x.t = 27rctd •^/[—7r,7r] for all 
combinations of x and t, where Ctd ^ [0,1] is a coefficient 
weighting dynamic disorder strength. To ensure that the QT 
distributions remain symmetric about the injection modes 
{x = 7,8) as both kinds of disorder are introduced, we mod¬ 
ulate 01 in the upper half of the array and 02 in the lower 
half of the array as shown in Fig. 3(a) (see Supplemental 
Information for details on transport symmetry). 

We perform a single QT experiment by initializing the 
PNP to implement a symmetric QW, programming a single 
realization (n) of a given strength of disorder, and mak¬ 
ing a measurement of the output distribution |^^(x,t/)|^ 
{tf = 7). Since a single realization of a level of disorder can 
result in a wide range of output distributions, this process 
is repeated over a large number of realizations {N). The 


FIG. 3. Quantum walk programs, (a) PNP phase settings for 
a QW with output-mode-symnnetric static and dynamic disorder 
response, (b) States are routed from time step ‘I’ to the output 
by programming all subsequent RBSs act as ‘wires’. 

ensemble average distribution is calculated as, 

1 ^ 

n=l 

We probe temporal evolution (measurements of (|^(x, t)|^) 
at each time step) by routing ‘wires’ (0 = tt, ^ = 0) for all 
RBS unit cells at time steps after the time step of interest 
as shown in Fig. 3(b). 

We will use the similarity metric^ 

where Dij {D-) is the experimental (computer simulated) 
distribution, to evaluate the closeness of our measured dis¬ 
tribution to a computer simulated distribution. 

Time-step tomography 

Figure 4 shows the measured spatiotemporal evolution of 
the quantum walker under four distinct environmental noise 
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FIG. 4. Spatiotemporal tomography, (a-d) (|^(x, t)|^) Vt (left) and (|^(x,t/)|^) (right), (a) Ballistic transport distribution 
with ctd = 0,ctid = 0. (b) Partially ballistic, partially localized distribution with ctd = ^^cud = 0.6. (c) Anderson localization with 
ctd = 0,cticz = 1. Grey line is a Laplacian fit to the distribution, (d) Incoherent, diffusive transport with ctd = ^^cud = 1- Grey line 
is a Gaussian fit to the distribution. 


conditions. We measure 200 unique unitary realizations for 
each level of disorder [ctd^cud) and time step (t^,l <i< 
7), totaling 5,600 experiments. 

First, we considered the case of an ordered QW 
(ctd, Qid)=(0, 0). In this case, we expect to observe a trans¬ 
port regime in which the distance of the quantum walker 
from the origin goes as 0{T), where T is the number of 
time steps taken. As shown in Fig. 4(a), ballistic transport 
is indeed what we observe, with similarity metric values for 
the mean distribution of St = 96.44% and Stf = 98.59% 
for all time steps and the final time step, respectively. Next, 
we consider the effect of an intermediate level of static 
disorder {ctd^ctid)={0, 0.6) on the output statistics. The 
measured distribution in Fig. 4(b) indicates intermediate 
levels of localization and ballistic transport characteristics 
{St = 94.86%, Stf = 99.40%). We observe that, by adding 
an intermediate level of static disorder, it is possible to gen¬ 
erate a more uniform distribution than would be possible 
for pure wave-like propagation or particle-like propagation. 
With strong static disorder (q^, Ctid)=(0, 1), the experi¬ 
mental data reveal an exponential localization to the initial 
sites (mode 7,8), as illustrated by the Laplace distribution fit 
in Fig. 4(c) with St = 96.85% and Stf = 99.68%. Finally, 
we consider the case of a fully incoherent quantum walk with 
Gid)=(l- !)■ The result is shown in Fig. 4(d) with a 
Gaussian distribution fit which highlights the particle-like 


transport statistics {St = 96.16% and Stf = 99.78%). 

To discriminate between the Laplace and Gaussian dis¬ 
tributions in Fig. 4(c,d), we computed the sum of squares 
residual errors as E = (where D is the mea¬ 

sured distribution and D' is the distribution fit) for both 
fitting distributions. A fit of (c) with the Laplace and Gaus¬ 
sian distribution gives E = 0.0028 and E = 0.0041, respec¬ 
tively: a fit of (d) with a Gaussian and Laplace distribution 
gives E = 0.0030 and E = 0.0057, respectively. In both 
cases, the distributions we chose were a better fit for the 
measured data. 


Environment-enhanced transport 

For the next set of experiments, we measure 400 lev¬ 
els of disorder with ctd.ctid ^ [0,1] and 200 unique QW 
realizations for each level, totaling 80,000 experiments. 
Since the QW distributions are symmetric about the in¬ 
jection mode {x = 8), we only need to consider modes 
1-7. For modes 1, 2, 13, and 14, (|^(x,ty)|^) is approx¬ 
imately constant across all disorder levels, as predicted by 
simulation, and so they are excluded from Fig. 5. Exper¬ 
imental measurements of (|^(x,t/)|^) for modes 3-7 (il¬ 
lustrated in Fig. 5(a)) are shown in Fig. 5(b-f). The 
measured data in Fig. 5(b-f) are in excellent agreement 
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FIG. 5. Complete transport space, (a) Rendering of output waveguides 1-14. We plot for modes 3, 4, 5, 6, 

and 7 which, due to symmetry, have the same behavior as modes 12, 11, 10, 9, and 8, respectively. We ignore modes 1, 2, 13, 
and 14 since they are constant across all ctd^cud parameters (due to the small number of possible paths that can constructively 
or destructively interfere on the edge of the light cone), (b-f) Experimental data corresponding to measurements of (|^(x,t/)|^) 
under 400 combinations of static and dynamic disorder. Fully coherent (incoherent), wave-like (particle-like) behavior occurs in the 
bottom left quadrant (right edge), (g) Transport efficiency {rf) through slice indicated by a dotted line along cud — 0.842 in (d.f). 
For a system with strong static disorder, the introduction of dynamic disorder enhances rj for outer modes 5 and 10 while reducing 
transport to modes 7 and 8. Each disorder level in the computer simulation corresponds to 20,000 realizations. 


with computer simulations, with final time-step similarities 
of Stf = 99.77,99.85,99.23,99.81,99.83% for modes 3-7, 
respectively. 

Ballistic transport, a signature of a wave-like random 
walk, occurs in ordered systems. As dynamic disorder is 
introduced, the quantum walker decoheres, losing its wave¬ 
like character: static disorder results in localization, increas¬ 
ing the transport efficiency to the injection mode. As shown 
in the lower left quadrant of Fig. 5 (b), this is indeed the 
general transport behavior at mode 3. The lower left quad¬ 
rant of Fig. 5(b) corresponds to ordered and nearly-ordered 
systems {ctd — O^cud — 0) in which a quantum walker 
starting at mode 7 or 8 is most likely detected at mode 3 
or 12. However, as disorder is introduced, the probability 
of detecting a quantum walker at x = 3 or 12 decreases. 
A similar behavior is observed for mode 4, as shown in Fig. 
5(c). 

The general behavior of modes 5, 6, and 7 shown in Fig. 
l(d-f) is opposite to that of modes 3 and 4 — increased 
static or dynamic disorder results in increased transport ef¬ 
ficiency to these modes with the exception of mode 5, which 
responds similarly for constant levels of dynamic disorder. 
In the regime of high static {cud > -6) and low dynamic 
{ctd < 0.4) disorder, the quantum walker is confined to 
modes 7 and 8 by Anderson localization, as indicated by 
the dark red region in the top left corner of Fig. 5(f). 

Environment-assisted transport occurs in (statically) dis¬ 
ordered systems in which transport has been arrested due 
to Anderson localization; through interactions with a fluc¬ 


tuating environment (dynamic disorder), the transport ef¬ 
ficiency of the quantum walker may be enhanced. Reben- 
trost et al.^^ predicted the existence of such an effect in 
continuous-time QWs. However, ENAQT has not yet seen 
experimental or theoretical consideration for DT QWs. 

Evidence of ENAQT for DT QW programs in the PNP 
can be found in Fig. 5(g). For strong static and weak 
dynamic disorder {cud — 0.8, Ctd — 0.), where Anderson 
localization is dominant, mode 5 exhibits a striking increase 
in transport efficiency with the introduction of dynamic dis¬ 
order. This enhancement in transport efficiency for mode 5 
corresponds to a reduction in the probability for the walker 
to stay at origin sites (x = 7,8), where the quantum walker 
is localized, as shown by the strong static disorder cut in Fig. 
5(g). In summary, a quantum walker confined to modes 7 
and 8 by Anderson localization is released through inter¬ 
action with a noisy environment, increasing the quantum 
walker’s transport efficiency to the outer modes. 

To determine whether a transport optimum that is simi¬ 
lar to the prediction for continuous-time QWs could be ob¬ 
served in a DT QW, we performed computer simulations of a 
PNP with 15 time steps. The data indicate that a transport 
optimum is indeed observable in systems such as the PNP 
(see Supplemental information for more information). 

DISCUSSION 

We have described the experimental observation of 
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environment-enhanced transport in discrete-time quantum 
walks over a set of 85,600 experiments spanning all combi¬ 
nations of static and dynamic disorder. This large number 
of realizations is enabled by a phase-stable, low-loss, large- 
scale, programmable nanophotonic processor implemented 
in a CMOS-compatible process. Each unit cell of the proces¬ 
sor can implement all SU(2) rotations in the spatial-mode 
basis; Reck et aP° explored a similar architecture for im¬ 
plementing the set of all single-particle unitary transforma¬ 
tions (U(N)). While single-particle quantum transport simu¬ 
lations are easily performed with classical computers, trans¬ 
port simulations exceeding several tens of photons^^ become 
intractable^^. Entangled photon sources^^’^^ and single pho¬ 
ton detectors^^ have recently been integrated into the silicon 
photonics platform, providing a path towards high photon 
number experiments. In addition, recently demonstrated 
low-latency superconducting logic devices^^ may provide 
fast on-chip feed-forward operations on quantum optical 
states. Together with programmable nanophotonic proces¬ 
sors, these chip-integrated technologies provide a promis¬ 
ing platform for future quantum simulation and computing 
tasks. 


METHODS 

Fabrication of the photonic circuit. We fabricated the PNP 
on a silicon-on-insulator wafer with 2 \im of buried Si 02 in the 
Optoelectronic Systenns Integration in Silicon foundry. Device 
regions are defined using partial- and full-etches in the silicon 
device layer resulting in 90 nm and 220 nm thick structures, re¬ 
spectively. A silica cladding 2\im thick is thernnally grown on top 
of the device layer to achieve transverse symnnetry for the optical 
modes. Input and output optical coupling is achieved with an in¬ 
verse taper from 500 nm wide to 200 nm wide over a distance of 
300p.m. Light is guided in 500 nm wide by 220 nm thick silicon 
ridge waveguides with a single transverse-electric mode (TEO) 
near 1550 nm, an effective index of 2.57^^. Thermo-optic phase 
shifters^^ are defined by full- and partial-etches and two Boron 
implants with concentrations 7-10^^ cm“^ and 1.7-10^° cm“^. 
We chose an electronic ‘acceptor’ dopant over a ‘donor’ dopant 
since the free-carrier absorption is less for free holes than free 
electrons"^^. Due to high active element densities, two aluminum 
layers and two via layers are required to simplify the phase shifter 
electrical signal routing. 
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